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FOREWORD 


This technical report is the final documentation of all work per- 
formed by Science Applications, Inc. for NASA/Marshall Space Flight 
Center under Contract No. NAS8-33022. The period of performance covers 
April -December 1978. Analysis methods and results presented herein are 
intended to assist NASA and DoE planners in assessing the viability of 
space disposal concepts within the overall context of nuclear waste 
management. 

The authors express their appreciation to the following who pro- 
vided valuable assistance in areas of technical discussion, data sources 
and report preparation: Gene Austin, NASA/MSFC Contract Monitor; Bob 

Nixon, NASA/MSFC; Jim Williams, OPL; Terri Ramlose and Kathy Osadnick, 
SAI Chicago Office. 
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LONG-TERM RISK ANALYSIS ASSOCIATED WITH 
NUCLEAR WASTE DISPOSAL IN SPACE 


1. INTRODUCTION AND SUMMARY 

Extraterrestrial disposal of hazardous waste material is being 
considered by NASA and DoE planners as a potential augmentation option 
of the national program for nuclear waste management. In support of 
this program planning effort, the present study addresses One aspect of 
the problem quite relevant to decision-making, namely, a quantitative 
measure of future risk and risk reduction requirements associated with 
the unplanned event of Earth reentry. Such an event, although rare, 
could result from either a failure occurrence in the deployment propul- 
sion system or long-term instability of a nominal storage orbit. 

The work reported herein is a continuation and logical extension 
of previous studies concerned with risk analysis [1,2]. These earlier 
investigations examined disposal destinations in heliocentric space and 
Earth-Moon space and were comprised of the following topical scope: 

(1) performance analysis of disposal options; (2) deployment accuracy 
and reliability; (3) orbit stability; (4) planetary encounter/Earth 
reentry probability; (5) long-term integrity of heat shield/containment 
materials; and (6) health consequences of unprotected reentry. 

With continuing in-house studies at MSFC [3], and their support con- 
tractors Battelle Columbus Labs [4] and Science Applications, Inc., a higher 
level of maturity has evolved in defining preferred disposal options and 
requirements. The current baseline concept would have waste payloads 
delivered to a stable solar orbit in the region between Venus and Earth; 
an alternate region of stability lies between Earth and Mars. Analytic, 
quantitative data that has previously been presented on questions of 
orbit stability and Earth reentry risk are thought to be basically cor- 
rect. However, as the space disposal concepts mature, it is important 
that this data be verified by more sophisticated (accurate) methods such 
as numerical integration and Monte Carlo statistical experiments. The 
present study then addresses this need for further verification analyses. 



1 .1 Study Objectives and Scope 

Three specific task areas of investigation were defined as described 
below. 

Task 1: Validation of Solar Orbit Stability . The objective of this 

task is to verify the quantitative bounds on the long-term (MO 6 years) 
stability of candidate storage orbits that were previously derived using 
secular perturbation theory. These orbits are initially circular and span 
the- range 0.83 to 0.88 AU between Venus and Earth, and 1.17 to 1.19 AU be- 
tween Earth and Mars. The verification analysis is concerned with evalu- 
ating the short-period perturbation effects which are ignored in the 
secular theory, and evaluating the analytical solution of the secular 
theory. This may be accomplished by numerical integration in the following 
manner: (1) short-term integration of the full equations of motion; and 

(2) long-term integration of the complete averaged equations of motion. 

Both methods are examined in terms of accuracy and computation time trade- 
offs. Numerical tests are performed for candidate storage orbits and sta- 
bility results compared with those of secular theory. 

Task 2: Validation of Earth Reentry Risk . The objective of this 

task is to verify the quantitative probabilities of planetary collision, 

and in particular Earth reentry, for "failed" payloads in solar orbit 

for time periods up to 5 * 10 6 years. Previous results were derived us- 
• • 

ing Opik's theory of planetary encounters to determine the orbit-conditional 
event probabilities which were then integrated over the deployment system 
failure distribution to obtain the overall risk profile. The new method- 
ology uses a Monte Carlo simulation of orbital evolution resulting from 
close planetary approaches to generate collision statistics as a function 
of time. Revised orbit-conditional event probabilities are obtained in 
this manner. The Monte Carlo procedure supersedes the "averaging" approach 
inherent to Opik's theory and presumably yields a more accurate statement 
of risk. A sufficiently large number of runs are made to generate statis- 
tical data for direct comparison with previous results. Revised risk 
profiles including retrieval (rescue) options are prepared as necessary. 


Task 3: Rescue Mission Requirements . The objective of this task is 

to evaluate the AV and propulsion stage requirements to carry out rescue 
missions (retrieval and placement) for a failed payload in solar orbit. 

The feasibility and requirements of automated rendezvous/docking in deep 
space are also assessed. 

It should be noted that the scope of the present study deals with 
the long-term aspects of risk assessment, i.e., failure events which occur 
after injection to deep space. It does not encompass systems failures and 
abort situations during the launch- to-orbit phase nor that portion of the 
early deployment phase while the payload is still in Earth orbit. These 
near-term failure modes are obviously important in the overall analysis 
of deployment safety and they are being treated in concurrent studies at 
MSFC and Battel le Columbus Labs. 

Results of this study are summarized briefly in the following sub- 
section. Organization of the main sections of the report follows the 
task sequence described above. Each section is independent for the most 
part and may be read separately in any order. 

1.2 Summary of Results 

Solar Orbit Stability . Numerical integration of the full equations 
of motion for time periods up to 12,000 years was performed for presumably 
stable orbits between Earth and Venus. The results indicate that secular 
perturbation theory is certainly valid (<10~ 3 AU error) for time scales of 
order lO 4 years and probably valid for much longer intervals. The sta- 
bility prediction of the analytic secular theory was then tested by numeric- 
ally integrating the complete averaged equations of motion for time periods 
up to 500,000 years. Figure 1-1 shows the orbit evolution comparison be- 
tween analytic and numerical secular theory. The agreement is very good 
in terms of the heliocentric distances covered and the general behavior of 
the orbital history. The increasing phase difference between the two cases 
is attributed to the truncation of higher order terms in formulation of the 




analytic theory. Similar results were found for an initially circular . 
orbit at 1.19 All. One interesting byproduct of the integration tests was 
the discovery of the possibly destabilizing effect of near-resonant orbits. 
The nominal 0.86 AU orbit is quite close to the 5:4 commensurability with 
Earth's orbital period, this occurring at a = 0.86177 AU. It would per- 
haps be best to select a new nominal target orbit that is slightly removed 
from this resonance, e.g., at 0.85 AU distance. Further work is recommended 
on the question of near-resonant orbit effects. As a general conclusion 
of the verification analysis, storage orbits between Venus and Earth or 
between Earth and Mars give every indication of meeting the desired sta- 
bility requirements for waste disposal, i.e., for time scales of the order 
10 6 years. The definitive test would, of course, involve numerical inte- 
gration of the complete equations cf motion for this period of time— a 
rather expensive research undertaking. However, it is suggested that this 
be done as soon as practical. 

Earth Reentry Risk . Analytic theory predictions of the long-term 
probability of Earth collision (resulting from deployment system failures) 
has been verified to within a close order-of-magnitude. A fortuitous re- 
sult is that the Monte Carlo data, presumably more accurate, indicate a 
lower risk by a factor of ^-5 for solar storage orbit missions. Figure 1-2 
shows a comparison of the analytic and numerical probability distributions 
for the nominal 0.86 AU destination orbit assuming a 99% reliability level 
for the deployment propulsion system. The revised results show that the 
reentry probability, for a single payload launch without rescue capability, 
is 6.4 x 10“ 4 for T £ 10 5 years, 2.2 x 10' 3 for T < 10 6 years, and 
5.7 x 10" 3 for T £ 10 7 years. Quite similar results are obtained for the 
solar orbit mission at 1.19 AU. Mission planners, quite understandably 
risk-adverse, have stated that a much smaller reentry probability of order 
10” 8 might be acceptable. This means, for example, that in a total space 
disposal program comprising 100 launches, there is only one chance in a 
million that any payload will reenter Earth's atmosphere within a time 
interval of order 10 6 years after launch. It appears that the only 
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COMPARISON OF ANALYTIC AND MONTE CARLO PREDICTION OF EARTH REENTRY RISK 
FOR SOLAR ORBIT MISSION (r = 0.86 AU, i « 2°) 




practical way to assure such low risk is through rescue mission capability. 
Figure 1-3 presents the revised 10 6 year risk profile for the 0.86 AU 
storage orbit; the effects of system reliability (R) and rescue mission 
degree of redundancy (N) are shown parametrically following the relation- 
ship P(R,N) = P(R,0)(1 - R)^. For example, with R = 0.99, each additional 
rescue mission yields 2 orders-of -magnitude further reduction in collision 
probability. The 10" 8 risk level can be attained with two rescue missions 
if the reliability is 0.996, or at most four rescue mission attempts if 
the reliability is only 0.965. 

Rescue Mission Requirements . A failed payload in solar orbit can be 
retrieved and placed into a stable orbit using the same propulsion hardware 
(OTV and kick stage*) employed for the nominal payload deployment. However, 
since the AV budget associated with rescue operations is always larger than 
nominal, it may be necessary in many instances to utilize two baseline kick 
stages. A possible alternative option is also available wherein a single 
kick stage is utilized to recover a stable orbit slightly different than 
the nominal target orbit. Table 1-1 presents one example from the range 
of failure cases studied. In all cases, the rescue mission is launched 
within one year of the nominal payload launch date (T^) , and stable orbit 
placement occurs within five years of that date. An important assumption 
apropos to the analysis is that the failed kick stage be jettisoned after 
rendezvous; it would be very inefficient to have to carry this dead weight 
through orbit transfer maneuvers. The most critical aspect of rescue mis- 
sions is undoubtedly the ability to rendezvous and dock with a payload in 
deep space. The preliminary assessment made here is that automated opera- 
tions are technically feasible based on current technology if a cooperative 
rendezvous mode can be assured . Such assurance implies a limited time 
interval between nominal payload launch and retrieval, and probably some 
degree of redundancy in payload vehicle systems (attitude control and 
communications) to enhance operational reliability. Some of the key fea- 
tures and requirements of rescue operations are summarized in Figure 1-4. 

*Kick stage (also referred to as SOIS - Solar Orbit Insertion Stage) is 
required for circularization burn at 0.86 AU on nominal mission. 
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NUMBER OF STANDBY 
RESCUE MISSIONS 



SYSTEM FAILURE LEVEL (1 - R) 

10 6 YEAR RISK PROFILE FOR SOLAR ORBIT DISPOSAL WITH RESCUE CAPABILITY, 
NOMINAL ORBIT: 0.86 AU CIRCULAR, 2° INCLINATION 




Table 1-1 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 


CASE NUMBER : 2B (Placement at T L + 4.240^ 

FAILURE EVENT : OTV and kick stage failure at T L 

Payload in orbit 0.93 x 1.0 AU, 1=2° 


RESCUE RESPONSE: 


Launch two* baseline kick stages for rendezvous, 
transfer and placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 

Time 

Distance 

C 3 or A V 

Inject to 0.906 x 1.0 AU 
orbi t 

Wait three revolutions 

t l + l.O 7 

1.0 AU 

C 3 = 1.61 0 1 = 2° 

Rendezvous and dock 

T l + 3.792 y 

1.0 

AV = 0.295 km/sec 

Transfer to 0.86 x 1.0 AU 

T l + 3. 792 y 

1.0 

0.598 

Placement 

T l + 4 . 240 * 

0.86 

1.187 

STAGE REQUIREMENTS: 

Stage Number 

A V (km/sec) 


Propellant (lbs) 

1 (Baseline Kick) 

0.893 


8829 (8397 used) 

2 (Baseline Kick) 

1.187 


8829 (7591 used) 

TOTAL INJECTED WEIGHT: 

2 x 11943 = 23886 

tl 

CO 

O 

O 

4-> 

to 

-Q 

1.61 


TOTAL RESCUE MISSION TIME: 3.240 years 


*Alternative using only single kick stage is placement into circular orbit 

0.888 x 0.888 AU, or elliptical orbit 0.86 x 0.917 AU. 
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Fig. 1-4 

AUTOMATED RENDEZVOUS AND DOCKING IN SOLAR ORBIT 
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2. VALIDATION OF SOLAR ORBIT STABILITY 


2.1 Background 

The extraterrestrial option for the storage of nuclear wastes requires 
orbits that are known to be stable, i.e., do not cross or come close to the 
orbit of any planet for a time interval around 500,000 years, which would 
allow long-lived radionuclides to decay through many halflives. Such an 
interval is very long when compared with orbit lifetime requirements for 
most manmade satellites, typically only a few years or decades. However, 
the 5 x lo 5 year interval is comparable to the longest period variation in 
the eccentricity of planetary and asteroid orbits while it is short com- 
pared with the orbital lifetime of planet-crossing asteroids and comets in 
the inner solar system (10 7 to 10 8 years). In an earlier study (Friedlander 
et al . [1]), the theory of secular perturbations for small bodies, origin- 
ally developed for investigations of the long-term motion of the asteroids, 
was applied to investigate the stability of various orbits initially located 
in the "gaps" between the terrestrial planets. The method of secular per- 
turbations requires that the short period perturbations of the orbit be 
eliminated by averaging over position in the perturbed orbit and the 
perturbing body's orbit. Analytic solutions for the orbit evolution can be 
obtained if the averaged equations are truncated after first order in the 
disturbing mass and second order in the eccentricities and inclinations 
(Brouwer and Clemence [2], Brouwer [3]). The resulting equations of motion 
reduce to a system of coupled, linear equations. This theory was used to 
determine the range of semimajor axes (a) between Earth and Venus and be- 
tween Earth and Mars, for which initially nearly circular orbits would be 
stable in the sense of avoiding neighboring planetary orbits for intervals 
of 500,000 years. The interplanetary gap between Venus-Earth and Earth-Mars 
together with the orbital evolution for representative orbits is given in 
Figures 2-1 and 2-2, taken from the earlier study. These initially cir- 
cular orbits at 0.86 and 1.19 AU with 2.0° initial inclination are adopted 
for now as nominal storage orbits. 


Earth Perihelion 



TIME (10 b YEARS) 

ORBITAL VARIATIONS OF AN INITIALLY CIRCULAR ORBIT AT a « 0.86 AU 




1.40 



TIME (10 6 YEARS) 

Fig. 2-2 ORBITAL VARIATIONS OF AN INITIALLY CIRCULAR ORBIT AT a = 1.19 AU 


The use of secular perturbation theory is necessary to survey many 
potential storage orbits over long intervals; however, there are various 
approximations made in the theory as well as restrictions on the applica- 
bility of the theory itself which give rise to questions regarding the 
accuracy of the analytic solution. Secular perturbation theory has been 
applied to studying the long-term motion of the planets and was given a 
partial vote of confidence by numerical integration. While the most 
recent analytic secular perturbation solution for the solar system (Brouwer 
and Van Woerkom [4]) includes all planets save Pluto, the numerically in- 
tegrated comparison by Cohen et al. [5] treated only the outer planets, 
Jupiter through Neptune. The inner planets were included only by the 
device of "casting them into the Sun," i.e., augmenting the Sun's mass by 
the mass of the terrestrial planets. In general, the agreement between 
secular theory and numerical integration was quite good with the principal 
discrepancies arising in the behavior of Uranus and Neptune. This discrep- 
ancy is probably due to Brouwer and Van Woerkam's ignoring terms in the 
analytic expansion which would account for the near 2:1 communsurability 
between Uranus and Neptune. Their secular theory neglects such resonance 
effects. No confirming numerical studies have been undertaken for terres- 
trial planets nor for application of secular perturbations to small bodies; 
however, there has not been any reason to doubt the secular results as long 
as resonances are avoided. Clearly, though, when orbits are being proposed 
for storage of nuclear wastes, it must be demonstrated that the stability 
analysis is valid. 

This study sought to confirm results from secular theory by numerical 
integration. Direct integration of the full equation of motion for a body 
in the inner solar system over a timescale of 10 5 to 10 6 years is an ex- 
pensive undertaking. A review of the literature of long-term integrations 
indicates that timescales of order 10 6 years have been covered, but only 
for the outer planets where the orbital periods are much longer. The 
integration of the outer planets to confirm secular theory covers 10 6 years. 



or 8.4 x 10 4 orbits of Jupiter. Other integrations such as the 4.5 x 10 6 
year integration of Pluto (Williams and Benson [6]) span only 1.4 x 10 4 
revolutions of Pluto. In the inner solar system, integrations cover 
roughly the same number of orbits, but give a shorter integration inter- 
val due to the shorter orbital periods, e.g., the orbit of asteroid 
1685 Toro was integrated for 5000 years (3125 orbital periods of Toro). 

An orbit between Earth and Venus with a = 0.86 AU has a period of nearly 
0.8 years; hence, a 500,000 year integration would require nearly 
6.3 x io 5 orbits, nearly an order of magnitude greater than has been done 
before. For this validation study, a two-part approach was adopted: 
first, to numerically integrate over a 500,000 year interval the complete 
averaged equations of motion to find the effect on the analytic solution 
of truncating higher order terms in the orbit elements. The second part 
consists'of numerically integrating the full equations of motion over as 
long an interval as possible given the available resources. 

2.2 Numerical Integration Program 

The integrations were performed using an integration program origin- 
ally developed by J. 6. Williams for studying the long-term motion of 
planets and asteroids. This program has the capability of integrating 
both the full equations of motion and the averaged equations of motion, 
including higher order terms neglected in the analytic solution. The pro- 
gram calculates the evolution of a non-resonant, non-intersecting orbit, 
which is perturbed by any of the planets Mercury through Neptune. Since 
we are interested in long-term behavior, the planetary orbits cannot be 
treated as fixed, but rather the orbital changes due to their mutual 
perturbations must be included. The Brouwer-Van Woerkom solution for the 
orbits of the eight principal planets is incorporated into the model. 

The program integrates the desired formulation for the variation of the 
orbit elements using a fourth order Runge-Kutta algorithm with either 
fixed stepsize or variable stepsize to keep the changes in any desired 
orbit element within specified bounds. The program presently is operational 



on the JPL Univac 1108 and performs integrations with double precision 
arithmetic, giving about 17 decimal digits of precision. The program takes 
about 0.1 CPU seconds per integration step of the full equations of motion 
when all perturbing planets are included. 

The long integration times dictate that compromises be made between 
integration accuracy and running time. Accuracy is improved by using 
smaller stepsizes, but smaller steps izes require longer execution times 
to .cover the same interval. Several runs were made over intervals of up 
to 2000 years for the full equations of motion and to 5000 years for the 
averaged equations to find the largest stepsize which does not incur a 
significant loss of accuracy. Table 2-1 compares the eccentricity from 
integrations of the full equations case with different stepsizes. Time- 
steps of the order of the orbital period or larger give results consider- 
ably different from short timestep runs, while there is little difference 
between the five-day and 0.05 year stepsize runs, and the 0.1 year and 
0.05 year cases are in reasonably good agreement. The significant loss 
of accuracy occurs for stepsizes greater than 0.1 years — clearly the 1 year 
stepsize is too large. Hence, 0.1 years seems the best stepsize for long- 
term integration of the full equations. Comparison of integration using 
this fixed stepsize with the variable step with Ae_ av = 3 * 10” 5 shows no 
significant loss of accuracy; the fixed stepsize case ran about twice as 
fast. 

For integration of the averaged equations of motion, the timestep 
can be considerably longer since we are no longer concerned with short 
period variations. The timestep can be longer than the orbital period, 
but should be small compared with the timescale of significant changes in 
the orbital elements. A series of 5000 year runs was undertaken to deter- 
mine the maximum timestep without significant accuracy loss; the results 
of this investigation are summarized in Table 2-2. Only stepsizes longer 
than 100 years are not acceptable while the 1 year step buys little addi- 
tional accuracy relative to a 10 year stepsize. The desired best stepsize 
is in the range of 10 to 100 years. Timing tests indicate a 500,000 year 



TABLE 2-1: ECCENTRICITY COMPARISON USING FULL EQUATIONS OF MOTION 

WITH DIFFERENT STEPSIZES 



Initial orbit: 

a = 0.86 

AU, e = 1 x 

10~ 5 , i = 

2.0° 

Time After 


Stepsize 

(years) 


Epoch {ye 

iars) 1 

.1 

variable 

i w 3xl0 ' D 

.05 

5 Days 

A) Long Term Integration 

0 lx 10 -5 1 x 10" 5 

l x io -5 

1 x IO -5 

1 x 10~ 5 

500 

.00474 

.00328 

- 

- 

- 

1000 

.01141 

.00410 

- 

- 

- 

1500 

.02150 

.00759 


- 

- 

2000 

.02625 

.00870 


- 

- 

B) Intermediate Term 





50 

- 

.00075 

.00058 

.00066 


100 

- 

.00127 

.00129 

- 


150 

- 

- 

- 

- 

- 

200 

- 

.00203 

.00204 

- 

- 

C) Short 

Term 





5 

- 

.000397 

- 

.000464 

.000472 

10 

- 

.000687 

- 

.000741 

.000748 

15 

- 

.00105 

- 

.00108 

.00109 

20 

- 

.00140 

- 

.00138 

- 

25 

- 

.00136 

- 

.00132 

- 

30 

- 

.00139 

- 

.00134 

- 

35 

- 

.00145 

- 

.00138 

- 

40 

- 

.00111 

- 

.00102 

- 

45 

- 

.00085 

- 

.00076 

- 

50 


.00075 

- 

.00066 
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TABLE 2-2 s ECCENTRICITY COMPARISON USING AVERAGED EQUATIONS WITH DIFFERENT 
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integration would take about 1.5 hours of computer time with the 100 year 
stepsize and proportionately longer with the smaller step. 

Further tests were made using the averaged equations of motion. 

Since the euqations integrate average elements, the initial elements must 
be average ones rather than osculating elements. Average elements were 
computed over a 20 year interval which is nearly the synodic period of 
Jupiter and Saturn. This should be the longest period significant per- 
turbation on non-resonant orbits in the inner solar system. The short 
period eccentricity variations are illustrated in Figure 2-3 for the 
calendar interval 1965 to 1995; these results were computed by integrating 
forward and backward 15 years from the 1980.0 epoch. The corresponding 
average elements are: 

a = 0.85993 AU 
e = 4.34 x Ur 4 
T = 1.999° 

The only significant change between the averaged and osculating ele- 
ments occurs in eccentricity where the average value is over 40 times the 
initial (osculating) value, but this is still essentially circular for 
applications purposes. As will be discussed in more detail later, there 
are indications that the nominal 0.8600 AU orbit is perhaps affected by a 
nearby resonance, namely, the 5:4 communsurability with Earth, which is 
exact at a = 0.861774 AU. To insure that the results were not modified 
due to proximity to the resonance, the nominal storage orbit was changed 
to one with a semimajor axis of 0.85 AU. An integration using this new 
orbit with initial elements e = 1 x 10“ 5 , i = 2.0° yielded mean elements 
of e = 1.46 x 10” 4 and T = 2.0003°; however, even in this orbit, the mean 
eccentricity is still much larger than in initial osculating value. A 
better initial orbit for comparison between integrated and averaged equa- 
tions is one that has initial osculating elements the same order of mag- 
nitude as the averaged elements. For semimajor axes around 0.85 AU, the 
mean eccentricity is typically (1-5) x 10" 4 . This value is comparable to 




1965 66 68 70 72 74 76 78 I960 82 84 86 88 90 92 941599 

FIGURE 2-3 ECCENTRICITY VARIATION FOR NEARLY CIRCULAR ORBIT AT 0.86 AU. 


the expected eccentricity resulting from guidance and execution errors 
associated with the insertion process which were shown to have a standard 
deviation of *\4 * 10“ 4 at this heliocentric distance (Friedlander et al . 
[1]). In light of these considerations, an initial osculating eccentricity 
of 5 x 1CT 4 was chosen which resulted in the average elements a” = 0.849975, 
e = 4.03 x 10 -4 ,- and i = 1.998. Integration accuracy tests were performed 
with the averaged equations using these initial conditions and are summar- 
ized in Table 2-3. The comparisons based on this orbit indicate little 
loss of accuracy in going from stepsizes of 10 years to 100-200 years. 

In light of these results, a 100 year timestep was used for numerical 
integration of the averaged equations. 

Numerical integration accuracies over longer intervals may be estimated 

based on analysis by Brouwer [7], who showed that numerical integration 

3/2 

errors accumulate proportional to n for the mean anomaly (position vari- 
able in the orbit), but for the other elements errors accumulate only as 
1/2 

n , where n is the number of integration steps. Assuming an error pro- 
1 /2 

pagation as n for the eccentricity, an eccentricity difference of 
1 x lo” 4 in 50 years (Table 2-1) would increase to 0.01 in 500,000 years. 

For circular orbits midway between Venus and Earth, an eccentricity of 
0.08 is required before the orbit can theoretically approach a planet. 

Hence, the integration error is apparently small compared to the eccentric- 
ity variations induced by planetary perturbations over this long time in- 
terval . 

2.3 Analysis 

Before commencing the long-term integration test runs, a 1000 year 
test run was made comparing the analytic averaged equations, integrated 
averaged equations and the fully integrated equations for the 0.86 AU ini- 
tial orbit. As seen from Figure 2-4, the analytic and integrated averaged 
equations have slightly different slopes, but are essentially in good 
agreement. This is true only over this limited interval. Clearly, if 
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TABLE 2-3: ECCENTRICITY COMPARISON OF AVERAGED EQUATIONS WITH DIFFERENT 
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this divergence continued over the full 500,000 years, there could be 
a major problem. However, as will be shown later, this divergence is due 
to a slight difference in the period of the secular variation but does not 
influence the amplitude. The full integration result plots eccentricity 
averaged over 20 years. For comparison, the results using the averaged 
equations, but starting with osculating elements, are given. 

The divergence of the integrated and averaged results on this short 
timescale is of some concern and two possible causes were explored: 

(1) proximity to the 5:4 resonance with Earth, or (2) semimajor axis vari- 
ations. A 0.85 AU orbit is well away from low order resonances with Earth 
or Venus. The 1000 year test run was therefore repeated using this orbit 
(Figure 2-5). The eccentricity variations for the fully integrated orbit 
are quite different than before, since the large amplitude, 60 year period 
oscillation has disappeared. The integrated elements now parallel the 
averaged elements quite well, with only a small divergence in the slopes 
and a small initial offset. In Figure 2-6, the semimajor axis variation 
of the 0.86 and 0.85 AU initial orbits are compared and show the much 
larger amplitude oscillations for the 0.86 AU orbit. These oscillations, 
which are likely introduced by the nearby resonance, have a mean period of 
nearly 60 years and probably are responsible for the eccentricity varia- 
tions with nearly the same period. The more stable 0.85 AU orbit clearly 
is a better choice for the validation study. 

A 500,000 year integration of the 0.85 AU orbit using the averaged 
equations is compared with analytic secular theory in Figure 2-7. In terms 
of the heliocentric distances covered and the general behavior of the 
orbital evolution, the agreement is quite good. The eccentricity ampli- 
tudes differ only slightly, and there is no indication that the integrated 
orbit will approach a planetary orbit in 500,000 years. However, a de- 
tailed comparison indicates a small difference in period of some terms in 
the eccentricity variation. The 25,000 year period variation is slightly 
longer in the integrated run than in the analytic evaluation; hence, there 
is an increasing phase difference between the two runs. These differences 
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TIME (10 years) OVER 500,000 YEARS 

FIGURE 2t 7 COMPARISON OF ORBIT VARIATIONS BASED ON ANALYTIC SECULAR THEORY (heavy 
lines) AND INTEGRATED AVERAGED EQUATIONS (thin lines) . 



are possibly due to the truncation of higher order terms in formulation 
of the analytic theory or to an error in one or both solutions. In gen- 
eral, however, there is no significant difference between analytic and 
integrated averaged equations of motion from the viewpoint of orbital 
stability. 

The most significant test of the averaged equations is to compare 
their results with those from integration of the full equations of motion. 
Figure 2-8 compares these results for the 0.85 AU orbit for a 12,000 year 
interval. Again, the overall agreement is quite good, indicating that 
secular perturbation theory does accurately describe the long-term evolu- 
tion of non-resonant orbits, at least on a timescale of 10 4 years. 

The period difference that was noted between the full averaged equa- 
tions and the analytic solution also appears in the fully integrated com- 
parison. The analytic solution shows a slightly shorter period than is 
indicated by the integrated solution, which is consistent with the trend 
noted earlier. This difference is of little significance for orbital 
stability. The important comparison, namely, the range of heliocentric 
distance covered due to planetary perturbations, shows excellent agreement. 

Comparison of integrated averaged euqations and analytic theory for 
an orbit between Earth and Mars is given in Figure 2-9. The initial orbit 
is circular at a = 1.19 AU with a 2.0° inclination. The orbital period is 
1.298 years which is between the 4:5 and 3:4 resonance with Earth, and 
based on the ratio of orbital periods with Mars (0.68969) is between the 
3:2 and 4:3 Mars commensurability. The exact 3:2 commensurabi 1 i ty with 
Mars (at a « 1.163 AU) is the closest of the two, so the 1.19 AU orbit is 
well away from any planet resonances. The comparison shown in Figure 2-9 
is generally good in that the heliocentric range is about the same between 
the two runs; however, the period difference still exists in the eccentri- 
city variations. An extended full integration would aid in deciding the 
source of this discrepancy. 
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2.4 Discussion 


The validation of secular perturbation theory in this study shows 
that the orbital extremes predicted by analytic secular theory are valid 
for timescales of order ICr years and probably for even longer intervals. 
There is no indication of higher order perturbations resulting from trun- 
cation of the full averaged equations in formulating the analytic solution, 
at least on timescales of order 10 6 years. Certainly from the perspective 
of waste disposal requirements, storage orbits between Venus and Earth or 
between Earth and Mars give every indication of meeting the desired sta- 
bility requirements. The definitive test would, of course, be to do a 
full numerical integration of the actual equations of motion for 
5 x 10 5 years. 

The above conclusion that the solar storage orbits do appear stable 
should be compared with results of other studies addressing similar ques- 
tions. One investigation of the stability of the solar system (Birn [8]) 
claims that there are no stable regions between Venus and Earth, while 
between Earth and Mars, the stable zone is from 1.16 AU to 1.39 AU. His 
method for stability analysis is to numerically integrate a four-planet 
system— Venus, Earth, Mars and Jupiter— for intervals of 700 to 1500 years. 
Birn employed the device of augmenting the planetary masses to shorten 
the timescale of dynamical interactions. This method assumes that if 
there exist secular variations of magnitude m J in an element (j, an integer, 
is the order of the perturbation), then increasing the mass of the per- 
turbing body (m) by a factor of f, will increase the magnitude of this 
term after time t, by a factor of f J . If an integration of T years is 
required to detect this secular variation due to m, then only an interval 
T/(f J ) is required for the augmented mass case. For example, if there are 
secular terms of order m 3 in the semimajor axis as discovered by Haretu 
in 1885*, an increase in perturbing mass of 10 would increase the magnitude 

♦More recent work by Message [9] indicates that this is just the leading 
term in a truly periodic function and that there are no known truly 
secular terms in semimajor axis to any order. 



of this secular term by 1000. If it takes a 10 6 year integration to de- 
tect secular terms using a perturbing mass m, they should be detectable 
in 1000 years with a perturbing mass of 10m. This method also assumes 
that other dynamical features of the problem such as the amplitude of 
periodic terms are not changed by the augmented mass. 

Birn found that orbits in the terrestrial planet gaps became unstable 
on accelerated timescales from 10 to 30 years when the masses of the ter- 
restrial planets were augmented by 100 and that of Jupiter by 10. However, 
there is no discussion of the behavior of individual orbit elements, par- 
ticularly a^ and e, nor of the physical nature of the instability. 

Presumably instability is caused by perturbations on the planetesimal 
orbit (a^ and/or e) causing it to become a planet-crossing orbit and to be 
eliminated by impact or ejection. It is well-known that there are no 
secular changes to a^ through second order in m but there are first order 
changes to e (hence the name secular perturbation theory!)- Increasing 
the masses of Earth and Venus by 100 would increase the magnitude of these 
terms by the same factor (assuming that Venus and Earth perturbations are 
dominant for an orbit midway between these planets). As was noted earlier 
and can be seen from Figure 2-1, eccentricities around 0.08 theoretically 
allow the planetesimal to become planet-crossing, while eccentricities 
twice as large guarantee that it will be planet-crossing. From Figure 2-5, 
the 0.85 AU orbit develops an eccentricity of 0.0008 on a timescale of 
300 real years and a value of 0.0016 in about 600 years. With a 100-fold 
increase in the perturbing mass, the eccentricity should be 0.08 to 0.16 
within 3 to 6 accelerated years, beyond which time a collision or close 
planetary encounter becomes a matter of random chance as described in 
Section 3. These times are reasonably consistent with Birn's range of 
10 to 30 accelerated years for instability. However, if the above sce- 
nario is the cause of instability, it indicates a flaw in the method for 
it is well-known that secular perturbation theory gives periodic variations 
in the eccentricity on long timescales (Figure 2-1). Increasing the initial. 



apparently secular change (Figure 2-5) to the point where the orbit be- 
comes planet-crossing changes the dynamical nature of the problem, i.e., 
the augmented perturbations seem to have introduced the instability. 

Note also that the first-order scaling would indicate instability on 
times of 1000 to 3000 years for the actual solar system whereas this study 
shows stability to 12,000 years for the present-day solar system. If the 
semimajor axis has secular variations, i.e., if there is a third order 
secular increase in a, then orbits should be stable for times of 10 4 to 10 7 
years (the fact that planetary masses were changed by different factors 
makes it somewhat difficult to scale the results). The lower number re- 
sults if Jupiter perturbations are dominant while the larger one holds if 
Venus and Earth are the largest perturbers. It should be pointed out, 
that for the outer planets, Birn argues that for a 100-fold increase in 
planetary masses, only the present-day semimajor axes of these planets are 
stable. However, in contrast to this conclusion, Nacozy [10] finds that 
Saturn would become unstable if the Jupi ter-Saturn system mass is increased 
by 29.25. It is clear that further work is required in this area to under- 
stand the physical nature of the all edged instability and the timescale on 
which it occurs. 

In summary, the present study shows by numerical integration that or- 
bits between Venus and Earth are definitely stable on timescales of 10 4 
years, while integration of the averaged equations shows that orbits be- 
tween Venus and Earth or between Earth and Mars are stable for at least 
5 x lo 5 years. No indications of incipient instability were found during 
the integration runs. Conflicting results regarding stability of orbits 
between Venus and Earth are attributed to differences in the dynamical 
models, and it is questioned whether results based on augmenting the masses 
of perturbing bodies are comparable to those presented in this study. 
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3. VALIDATION OF EARTH REENTRY RISK 


The field of statistical celestial mechanics, developed in the 1950's 
by the pioneering work of Ernst Opik, is now recognized as an important 
contribution to scientific inquiry of solar system evolution and pro- 
cesses. Its principal application has been to studies of asteroid colli- 
sions, planet cratering and related cosmogonic processes. The need for a 
statistical representation of actual events (past or future) derives from 
the practical limitations of classical celestial mechanics to definitively 
model orbital evolution over millions of years. This is especially the 
case when such evolution involves strong gravitational interaction result- 
ing from repeated close planetary encounters. Opik's analytic theory, 
based on geometrical reasoning and convincing arguments of averaging, is 
simple in form and allows rapid evaluation of event statistics. 

A nuclear waste payload which fails to achieve its "stable" destina- 
tion orbit may be considered a stray body in the solar system. Stray 
bodies in interplanetary space which cross and possibly intersect the 
orbits of planets define a long-term process of probable orbital change 
and eventual elimination through the events of planet collision, solar 
impact or solar system escape. The problem of quantifying the probabil- 
istic fate of a failed payload is identical in all respects to the general 
problem of statistical celestial mechanics. Opik's theory was utilized 
in the development of the PEPA computer program from which previous study 
results were generated [1]. These results are thought to represent 
"statistical truth" to within a close order-of -magnitude, i.e., a factor 
of 2 to 5. It is the verification of this presumption which is the sub- 
ject of this section of the report. 

3.1 Monte Carlo Simulation of Planetary Encounters 

Since Opik's formulae are based on averaging approximations to real 
world orbital perturbations, certain limitations of the analytical approach 
arise when dealing with complex problems. The main difficulty is related 
to the following two factors: (1) the dispersion in relative encounter 



velocity magnitude caused by successive close encounters with a real 
planet moving in a precessing, elliptical orbit rather than a fixed cir- 
cular orbit; and (2) the dispersion caused by multiple close encounters 
with more than one planet resulting in a "playing ball with the particle" 
effect. This dispersion or so-called acceleration of encounter speed is 
not accounted for in the basic theory. Although Opik in his later work 
recognized and explained the acceleration effect, after it was discovered 
empirically, it is difficult to ammend the formulae in a simple, straight- 
forward way to properly account for this effect. In contrast, the problem 
is readily treated by numerical simulation using Monte Carlo statistical 
methods. 

Other science researchers, most notably Arnold [2,3] and Wetherill 
[4,5], have extended the treatment of statistical celestial mechanics by 
adopting the numerical approach and we follow their lead here. The re- 
suit is a semi analytic method in that Opik's equations, which give the 
probability of close planetary encounter, are used in an iterative manner 
to drive a Monte Carlo generator of orbital evolution and collisional 
statistics. The method proceeds as follows: 

1. Input initial orbit elements (a,e,i) at time T. 

2. Identify all planets whose orbits may be crossed 
by the stray body. 

3. Using Opik's equations, calculate the probability 
of close approach (within K collision radii) to 
each planet. A typical choice is K = 5. 

4. Combine these probabilities to calculate a mean 
time (t) between encounters. 

5. Choose the next planet to be encountered, at 
random, weighted according to the relative en-. 
counter probabilities. 

6. Choose a time interval AT to the next encounter. 

at random, from an exponential distribution e" AT /x. 

Record the running time T T + AT. 


7. Choose a random location on the encountered planet's 
orbit within the allowable region of orbit crossing. 

8. Choose a random point on a target disk of K collision 
radii centered on the chosen planet. The stray body 
is assumed to move on a hyperbolic relative trajectory- 
passing through this point. The region (K/2,K) is 
weighted more heavily than the interior region 
(0,K/2) so as to account for those passages outside 

of K radii which may occur but are ignored for reasons 
of computational efficiency. 

9. If collision does not occur, calculate the change in 
orbit elements produced by the close gravitational 
encounter. 

10. Use the new orbit elements to reinitialize Step (1). 

11. Continue this procedure until some event terminates 
the evolution process (planet collision, solar impact 
or solar system escape). 

12. Upon termination of a given Monte Carlo pass, recycle 
to the original orbit state and repeat the entire 
random experiment. Statistics are accumulated over 

a large number (several hundred) of passes. 

This methodology may be recognized as a statistical analog of a multi - 
revolution numerical integrator. It also involves approximations to 
reality but at a higher level of accuracy compared to the strictly 
analytic, closed-form theory. Also, the results could be sensitive to 
the sequence of random number generation and the total number of Monte 
Carlo passes. As in any finite sampling experiment, care must be taken 
to generate a sufficient number of cases so as to have reasonable confi- 
dence in the statistical data. In practice, this implies the usual trade 
off between statistical variability and computational cost. 

The Monte Carlo simulation outlined above was programmed as a modifi 
cation to the existing computer program PEPA. The new subroutine may be 
selected upon user option in place of the analytic formulae for the com- 
putation of event probabilities conditional upon the specified initial 
orbit. A variety of test cases were run with reference to data published 



by Wetherill to assure that the program worked properly. Application 
was then made to the waste disposal problem as described in the remainder 
of this section. 


3.2 Comparison of Conditional Probability Results 


Before presenting example results, it will be helpful to review the 
basic definitions of "conditional probability" and "total probability." 
Conditional probability of an event E refers to the likelihood that the 
event will occur within a specified (long) time interval T given that the 
payload is left in an unstable orbit having initial conditions (a,e,i;t) 
resulting from a deployment system failure at time t. The point of 
failure occurs somewhere in the (short) time interval (0,t^) comprising 
the nominal deployment sequence. Total probability then refers to the 
likelihood of this same event as obtained by integrating the entire spec- 
trum of conditional probability over the time distribution of deployment 
system failures. Let R(t) denote the reliability function for the deploy- 
ment system. Total probability may be computed from the integral expres- 
sion 


p e (t) 


c f 

l p ‘ 


(T|t) 


dR(t) 

dt 


dt 


where P£(T|t) is the conditional probability of the event. Staged pro- 
pulsion sequences and startup reliability are readily accommodated in the 
definition of R(t). The PEPA program uses Simpson's Rule integration to 
evaluate total probability. 

In assessing the validity of previously obtained analytic results, 
it is the conditional probability P £ (T | t) that is of chief concern and, 
in particular, the event of Earth collision or reentry. In other words, 
does Monte Carlo simulation give essentially the same results for Pr(T|t) 
as predicted by Opik's theory, or is there a significant discrepancy? 

We have reason to expect some difference because of the encounter speed 
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dispersion effect that was mentioned earlier. This dispersion is par- 
ticularly important when the initial value is small, e.g., less than 0.1 
of Earth mean orbital speed. This is in fact the situation of interest 
for an Earth-crossing orbit having dimensions close to that of the Earth's 
orbit and inclined to the Earth’s orbit by a small angle. 

Consider an initial orbit of size 0.86 * 1.0 AU with a 2° inclination 
to the ecliptic plane. The Earth is the only planet that is crossed 
initially . However, the first velocity dispersion factor opens an effi- 
cient loophole to the eventual occurrence of Venus crossing. When that 
happens the second dispersion factor comes into play with eventual cross- 
ing of Mars, then Jupiter. Once Jupiter is crossed its enormous gravita- 
tional field, on close encounter, can easily cause the object to be ejected 
from the solar system. Hence, we are dealing with a dynamic process that 
can lead to many different kinds of elimination events even though Earth 
collision may be the only event that can occur initially. This process 
of velocity dispersion or acceleration takes time to build up through 
numerous close but non-colliding encounters with Earth. The key question 
then from a probabilistic point-of-view is whether the time constant of 
this dynamic process is short or long compared to the time interval 
(e.g., T <_ 10 6 years) of concern for waste disposal risk. If the time 
constant is long, then we can expect fairly close agreement with the 
analytic prediction. 

A computer printout of results is shown in Table 3-1 for the condi- 
tional orbit resulting from failure of the circularizing burn at 0.86 AU. 
Opik's theory predicts a mean lifetime against Earth collision of 
t = 0.52 x 10 6 years with the probability being exponentially distributed 
as 1 - exp(-T/r). The Monte Carlo simulation shows widely different re- 
sults for elimination lifetime. Earth collision remains the most likely 
event of elimination, but only 55.8% of all cases result in this disposi- 
tion. Furthermore, the mean time of Earth collision is 16.4 x 10 6 years 
which is substantially longer than the analytic value. Venus collision 
is the next most likely event at 29.2% followed by solar system ejection 
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(due to Jupiter and Saturn crossings) at 10.6%. The combined mean life- 
time against all possible events of elimination is 37.3 * 10 6 years. 

Note that the mean value of Earth encounter speed at the time of collision 
is 0.172 EMOS which indicates the substantial dispersion from the initial 
value of 0.062 EMOS. 

The first conclusion we may draw from this test is that the disper- 
sion effect, not accounted for by the analytic theory, acts to increase 
the lifetime for this particular initial orbit. Although this result is 
certainly desirable in that the probability of Earth collision is de- 
creased, the real significance depends on the time interval of interest. 

The time profile is shown in Figure 3-1 which compares the Monte Carlo 
statistics with the analytic value of Earth collision probability. The 
agreement is fairly close for T < 10 5 years. The largest difference, at 
T -v 10 6 years, is only a factor of 3.8. Hence, the statement that analytic 
theory should predict the Earth collision probability to within a close 
order-of -magnitude appears to be correct. 

An analysis of Monte Carlo sampling variations is presented in 
Table 3-2 for the conditional orbit 0.86 x 1.0 AU. Earth collision prob- 
ability values are listed for a range of time intervals showing statis- 
tical data variations over five sampling cases each for 100 and 500 samples 
per case. As expected, the data accuracy (standard deviation) improves 
with larger sample size and longer time interval. With 100 samples, the 
standard deviations about the mean value are 44.8%, 25.7% and 13.6%, 
respectively, for time intervals of 10 5 , 10 6 and 10 7 years. Corresponding 
results for 500 samples are 13.3%, 5.1% and 4.7%. These variations are 
sufficiently small to place reasonable confidence in the results of a 
finite sampling experiment, i.e., as it pertains to characteristic dif- 
ferences between the analytic and numerical methods. 

A comparison of the two methods was also made for a very different 
initial orbit, namely, one with perihelion distance near Earth and aphelion 
distance near Jupiter. In this case the dominant event of elimination is 
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ejection out of the solar system due mainly to Jupiter encounters. The 
probability of Jupiter collision is about an order of magnitude less than 
ejection but much greater than any other planet collision event. Results 
are shown in Figure 3-2 where the probability of occurrence of the two 
dominant events are plotted as a function of time over the interval 
10 5 to 5 x lo 7 years. Monte Carlo and analytic results are in very good 
agreement for ejection and somewhat less so for Jupiter collision — but still 
within a factor of 3. Even though Earth's orbit is crossed by this ini- 
tial orbit the probability of Earth collision is relatively low. Analytic 
theory predicts a value of 10“ 4 for T = 5 x 10 6 years. In 1500 Monte Carlo 
samples only one Earth collision event was observed. Although a much 
larger sample size would be needed in this case for an accurate comparison 
with analytic theory, the limited data obtained does not indicate any 
major discrepancy. 

3.3 Comparison of Total Probability Results 

The following three waste disposal mission examples were selected for 
purposes of risk verification: 

1. Solar storage orbit, 0.86 AU, 2° inclination 

2. Solar storage orbit, 1.19 AU, 2° inclination 

3. Solar system escape, 0° inclination. 

The propulsion system assumed for the two cases of solar orbit deployment 
is a reusable OTV at injection and a kick stage for circularization at 
the nominal heliocentric distance (see Section 4). Although the solar 
system escape mission is not currently considered a baseline option, it 
is included here for completeness since it had been exampled previously [1]. 
A three-stage (large solids) propulsion system is assumed for this case. 

The propulsion system failure model, taken to be identical for each 
stage, includes a startup component (reliability R^j) and an operational 
component (end-point reliability R Ep ). For purposes of the risk verifi- 
cation analysis we will simply assume that R ST = R Ep = 0.99. The effect 
of the system reliability level will be described parametrically later in 
Section 3.4. 




We examine first the case of the 0.86 AU solar orbit destination. 

The range of conditional failure orbits encompasses perihelion distances 
from 1.0 to 0.86 AU and aphelion distances from 1.0 to 0.86 AU. Total or 
integrated probability data are presented in Table 3-3 which compares 
analytic theory and Monte Carlo computer printouts. Earth collision is 
the only event of elimination predicted by analytic theory with probabil- 
ity varying from 2.7 x 10" 3 for T = 10 5 years to an asymptotic upper limit 
of 1.3 x 10" 2 . Monte Carlo results verify that Earth collision is the 
dominant event for shorter time intervals, but that other events, notably 
Venus collision and solar system ejection, become increasingly as likely 
as time grows. Figure 3-3 compares the probability versus time distribu- 
tion for Earth collision— the event of primary interest to disposal risk 
assessment. Monte Carlo results show a probability variation of 
6.4 x 10" 4 to 5.7 x 10" 3 over the time interval 10 5 to 10 7 years. This 
represents a reduced risk relative to analytic theory prediction by a 
factor of 5 or less. These results are consistent with the conditional 
probability comparison discussed earlier and with the accuracy bounds 
expected of the analytic theory. 

Similar data are presented in Tables 3-4 and 3-5 and Figures 3-4 
and 3-5, respectively, for the 1.19 AU solar orbit and the solar system 
escape missions. The probability characteristic and comparison with 
analytic theory of the 1.19 AU case is nearly the same as for the 0.86 AU 
disposal destination. Monte Carlo results indicate a reduced risk of 
Earth collision by a factor of 5. For the solar system escape destina- 
tion, a much closer agreement is found between analytic and numerical 
methods. Earth collision probability predicted by analytic theory varies 
from 2.5 x 10" 4 to 7.5 x 10" 3 over the time interval 10 5 to 10 7 years. 
Monte Carlo results are 80% higher at T = 10 5 years, approximately the 
same at T = 10 6 years, and 40% lower at T = 10 7 years. 



3.4 Risk Reduction with Rescue Mission Capability 

It seems appropriate to revise the waste disposal risk profiles based 
on the new data obtained from the Monte Carlo analysis. We shall examine 
the parametric effects of deployment system reliability and risk reduction 
through retrieval of failed payloads. With reference to the 0.86 All solar 
orbit mission, we have found that the probability of Earth collision is 
^2 x 10~ 3 over a time interval of 10 6 years. This means in effect that 
one payload in 500 launches can be expected to return and reenter Earth's 
atmosphere sometime within 10 6 years after launch. Let us suppose for the 
sake of argument that a decision maker is risk-adverse and therefore con- 
siders this level of risk to be too high. What means are available to 
reduce the risk? Improving reliability is the most direct way since de- 
ployment system failures are responsible for the problem in the first place. 
System reliability enhancement by design or redundancy is certainly pos- 
sible, but there is some practical limitation to this approach beyond 
which it is no longer technically or economically feasible. We have al- 
ready assumed a nominal reliability level of 99%. Increasing this to 
99.9% will only effect a ten-fold reduction in risk. It seems clear that 
a different approach is needed to obtain substantial improvement — say to 
a probability level of 10“ 6 or better. Payload retrieval or rescue capa- 
bility can provide this assurance. 

Retrieval is defined here as the ability to send another propulsion 
system to rendezvous with the failed payload in solar orbit and to place 
it into the desired stable orbit. Once this capability exists we may as 
well adopt a policy of redundancy on a mission level. That is, if the 
first rescue mission fails we launch a second mission, if the second fails 
we launch a third, etc. Therefore, assume a standby, multiple rescue 
capability, each having an equal probability of success or reliability R. 
The term "standby" does not necessarily imply an on-the-pad readiness since 
the rescue missions may be carried out over a period of several years. 

The chance of reentry during this early time period is essentially zero. 
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Fig. 3-3 COMPARISON OF ANALYTIC AND MONTE CARLO PREDICTION OF EARTH REENTRY RISK 
FOR SOLAR ORBIT MISSION (r = 0,86 AU, i = 2°) 
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Fig. 3-4 COMPARISON OF ANALYTIC AND MONTE CARLO PREDICTION OF EARTH REENTRY RISK 
FOR SOLAR ORBIT MISSION (r - 1.19 AU, 1 = 2°) 
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SOLAR SYSTEM EJECTION 




TIME INTERVAL T (YEARS) 

Fiq 3-5 COMPARISON OF ANALYTIC AND MONTE CARLO PREDICTION OF EARTH REENTRY RISK 
FOR SOLAR SYSTEM ESCAPE MISSION (i = 0°) 



Figure 3-6 shows the 10 6 year risk profile as a function of system 
failure level (1 - R) and the number of "standby" rescue missions (N). 

The data are calculated by the formula 

P ec (R,N) = P ec (R,0)[1 - R] N 

where P^RjO) is the total probability of Earth collision for a single 
payload launch obtained from the Monte Carlo analysis. For example, with 
R = 0.99, each additional rescue mission yields 2 orders of magnitude 
further reduction in collision probability. The shaded area in Figure 3-6 
indicates the region of likely interest for mission design. Suppose that 
the maximum acceptable risk of Earth reentry is 10~ 8 for a single payload.* 
This risk can be attained with two rescue missions if the reliability 
level is 0.996, or at most four rescue missions if the reliability is only 
0.965. Figure 3-7 shows the risk profile as a function of time over the 
interval 10 5 to 10 7 years assuming a system reliability level of 99 %. 

In summary, it has been shown that it is not necessary to have deploy- 
ment systems of extremely high reliability in order to achieve very low 
risk. Instead, rescue mission provides this assurance through the powerful 
mechanism of redundancy. It is suggested that such redundancy on a mission 
level is likely to be far less costly to obtain than the equivalent redun- 
dancy that would be required on a system design level. 


*e.g., in a total space disposal program consisting of 100 launches, this 
long-term risk corresponds to only one chance in a million that any pay- 
load will reenter Earth's atmosphere. 
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4. RESCUE MISSION REQUIREMENTS 

For the case of waste disposal in solar storage orbits, the proba- 
bility that a single failed payload will reentry Earth's atmosphere 
within a time period 250,000 years after launch is of the order 10" 3 , 
i.e., one chance in a thousand. This relatively high level of risk has 
been stated to be unacceptable by some program planners. However, the 
analysis described in the previous section has shown that rescue mission 
capability may be used to reduce this risk to less than 10 -8 . The ques- 
tion to be addressed then relates to the requirements imposed on such 
rescue missions. 

Consider the following situation and response scenario: 

• Situation- -due to some deployment system failure, the payload 
is stranded in an Earth-crossing orbit having a finite proba- 
bility of reentry in the long term. The payload has not broken 
apart and it appears likely that a rendezvous/docking operation 
can be achieved in the short term, i.e., within several years 
after launch. 

• Response— as soon as is practical, launch a rescue mission to 
retrieve and place the failed payload into a stable solar or- 
bit. If for some reason the rescue mission fails, then launch 
another— (etc. ) , thereby assuring the acceptable level of risk 
against reentry. 

Problem areas for study implied by this scenario are: (1) the range of 

possible failed orbits; (2) rendezvous phasing orbit requirements with 
rescue time/AV trades; (3) rescue propulsion stage requirements with 
possible constraints on hardware commonality; and (4) automated rendez- 
vous/docking procedures and guidance accuracy. 

Figure 4-1 serves to scope the analysis by describing the failure 
events considered and the range of resulting solar orbits. The desired 
target orbit is circular at 0.86 AU and inclined to the ecliptic plane 



by a small angle, nominally 2°. Injection to the 0.86 x 1.0 AU transfer 
orbit is performed by a reusable OTV propulsion system, and circulariza- 
tion at 0.86 AU (6 months later) is accomplished by a kick stage injected 
with the payload. Three failure modes occurring after Earth escape con- 
dition* (C 3 = 0) can leave, the payload in an unstable solar orbit. In 
the first instance, a late OTV failure, the proper response is to imme- 
diately start up the kick stage to implement completion of the nominal 
injection to the 0.86 x l.o AU transfer orbit. For corrective maneuvers 
within 20,000 km altitude above the Earth, the required aV is less than 
100 m/sec. Since this is within the propellant budget of the baseline 
kick stage, there is no further reason to consider this particular failure/ 
corrective mode— assuming of course that the kick stage does not fail. 

The two remaining cases of interest then are: (1) sequential failure of 

both OTV and kick stages at injection leaving the payload in an orbit 
with aphelion distance of 1.0 AU and perihelion distance somewhere in the 
range 1.0 to 0.86 AU depending on the time of failure; and (2) successful 
OTV burn but kick stage failure during the circularization maneuver leav- 
ing the payload in an orbit with perihelion 0.86 AU and aphelion somewhere 
in the range 1.0 to 0.86 AU. 

In the first case above, a rescue mission is clearly needed since 
both nominal propulsion stages are inoperable. In the second case, a 
rescue mission is required if the kick stage fails at startup, but may 
not be needed if the aphelion distance is sufficiently close to 0.86 AU 
as a result of a late kick stage failure. Rather than restrict the 
analysis by such implied operational decisions, we will examine the res- 
cue aV requirements over the entire range of orbit dimensions. However, 
we will return to the question of orbit stability for situations of off- 
nominal, but non-Earth crossing solar orbits. 


*Failure occurrences prior to achieving Earth escape will leave the waste 
payload in an Earth-bound orbit. Rescue requirements for this situation 
are being studied by MSFC and Battel le. 


Fig. 4-1 

FAILURE EVENTS AND RESULTING ORBITS 
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4.1 Rendezvous Phasing Orbits 


The first step in the analysis is to calculate the minimum AV 
requirement assuming a time-independent optimal phasing situation and 
no additional loss due to finite thrust penalty or terminal rendezvous 
maneuvers. This establishes a reference data base against which real 
maneuver requirements can be assessed. Results are shown in Figure 4-2. 

The first velocity impulse, &V ls is provided by the rescue OTV at Earth 
parking orbit injection. Over the range of orbits at 2° inclination, 
the AVj requirement is approximately constant (3.25 to 3.31 km/sec); 
the equivalent launch energy C3 range is 1.08 to 2.35 km 2 /sec 2 . The 
remaining impulses, aV 2 and AV 3 , are provided by the rescue kick stage. 

The right-hand graph which depicts failure mode #1 shows an ideal "free" 
rendezvous and a placement AV varying from the nominal value of 
1.187 km/sec to 2.330 km/sec. Failure mode #2 is shown in the left-hand 
graph. In this case, the sum aV 2 + aV 3 is always equal to 1.187 km/sec; 
i.e., the total AV requirement is independent of aphelion distance. To 
summarize, over the range of possible rescue orbits the minimum total AV 
requirement varies from 4.5 to 5.7 km/sec depending on the orbit size. 

It is possible to approach the ideal phasing conditions as closely 
as desired if one is willing to wait the requisite amount of time before 
mounting the rescue mission. Although it is true that the near-term risk 
of Earth reentry is negligible, there may be a practical desire to mount 
an early rescue mission; e.g., to take advantage of a still-functioning 
communications and attitude control system on the failed payload in order 
to cooperate in the rescue operations. It is important, at least, to 
identify the AV tradeoff with time. Figure 4-3 shows such a tradeoff for 
a failed payload in a 0.86 * 1.0 AU orbit with arbitrary inclination. 

Results are given in terms of AV penalty above the minimum AV requirement, 
but again assuming no losses other than orbit phasing. The left-hand 
graph applies to rendezvous and placement at 0.86 AU, while the right-hand 
graph assumes rendezvous at 1.0 AU and placement M).45 years later at 
0.86 AU. Data is shown for launching the rescue mission from 0.5 to 2 years 
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Fig. 4-2 MINIMUM (OPTIMAL PHASING) AV REQUIREMENTS FOR SOLAR ORBIT RESCUE 
AV, = PARKING ORBIT INJECTION; AV 7 = RENDEZVOUS; AV, = PLACEMENT 






(oas/ui>|) wnWINIW 3A09V AilVNBd AV 



RESCUE PLACEMENT (YEARS AFTER Tj RESCUE PEAEEMENT < VEARS AFTER V 

Fig 4-3 REAL PHASINS 4 V PENALTY FOR SOLAR ORBIT RESCUE (T L = PAYLOAD LAUNCH DATE) 

ASSUMING PAYLOAD ORBIT = 0.86 * 1.0 AU, PLACEMENT ORBIT = 0.86 AU CIRCULAR 




after the original payload launch date, and placement time up to 5 years 
after payload launch. It should be noted that rescue launch at T L + 0.5 y 
allows for only a 3 week period between first indication of kick stage 
failure and rescue launch, implying an "on-the-pad" readiness. This 
short response is not inconsistent with the possible need to quickly 
rescue a failed payload in Earth orbit. Another point to note is that 
only the solid dot data points correspond to specific rescue mission 
scenarios. 

Launches at half-year intervals are more efficient if rendezvous 
takes place at 0.86 AU, whereas a 1.0 AU rendezvous is better for launches 
at 1 year intervals. Also, for placement within 5 years, it is best to 
launch the rescue mission early, i.e., within the first year. The ear- 
liest placement, at 1.35 years, requires an immediate rescue mission hav- 
ing a aV penalty of 2.402 km/sec. Otherwise, the rescue mission launched 
at T l + l y with rendezvous at 1.0 AU provides a minimum A V penalty from 
1.693 to 0.369 km/sec. 

4.2 Propulsion Stage Requirements 

Rescue mission AV requirements are translated to propulsion system 
performance using the mass characteristics data provided by MSFC. This 
data listed in Table 4-1 assumes that the rescue mission is implemented 
by propulsion hardware similar to that used for the nominal payload 
launch. In most cases, a reusable 0TV will be employed for the rescue 
injection and one or more kick stages used for implementing the rendez- 
vous, orbit transfer and placement maneuvers. The nuclear waste payload 
is assumed to weigh 10,495 lb. Performance calculations include a 
100 m/sec allowance for terminal rendezvous maneuvers and finite thrust 
losses, and assume that use of attitude control propellant is equally 
split between pre- rendezvous and pre-placement maneuvers. Another impor- 
tant assumption is that the failed kick stage is jettisoned after rendez- 
vous; it would be very inefficient to have to carry this "dead weight." 



Table 4-1 

MASS CHARACTERISTICS DATA 
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Several failure scenarios covering the spectrum of Earth-crossing 
orbit conditions were postulated and details of possible rescue responses 
were calculated in terms of the aV and propulsion stage requirements. 
Results of these failure case studies are presented in Table 4-2. Rescue 
stage requirements shown reflect the desired placement into the nominal 
circular stable orbit at 0.86 AU within 5 years of the original payload 
launch. This may necessitate the use of two baseline kick stages as 
noted in Cases 1 (Alternative), 2A, 2B and 3A. It is always possible, of 
course, to design a single, larger propulsion stage to accomplish the res- 
cue mission. However, this implies the use of nonstandard hardware and 
may be less cost-effective than using two standard kick stages. 

One exception to using the standard reusable OTV is shown for failure 
Case 1 in which the payload is left in a 1.0 x 1.0 AU orbit. In this case 
an expendable OTV is considered for rendezvous and transfer AV implemen- 
tation, and a single kick stage performs the final placement maneuver. 

The cryogenic propellant in the OTV must be stored for a longer than us- 
ual time interval, but probably not exceeding 10 days of operation in the 
near vicinity of Earth. 

There is an attractive alternative to consider when a single base- 
line kick stage is insufficient to recapture the nominal 0.86 AU circular 
orbit. That is, the available AV capability could be used to establish a 
circular orbit at some other distance or a slightly elliptical orbit. 

Data for this alternative option are presented as footnotes in Table 4-2. 
If these fallback orbits exhibit a sufficient measure of stability against 
planet perturbations, then it should be possible to perform all rescue 
missions with a single baseline kick stage. Figure 4-4 shows the varia- 
tion of orbit distance for an initially circular orbit at 0.90 AU which 
encompasses the fallback options noted for Cases 2A and 2B. The storage 
orbit aphelion distance does barely cross Earth's perihelion, but not 
until a time interval of 600,000 years has elapsed. Hence, the 0.90 AU 
circular orbit could be said to possess at least marginal long-term 
stability. Additional results are shown in Figures 4-5 through 4-7 



Table 4-2 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 


CASE NUMBER: 


.1 


FAILURE EVENT: 


OTV and kick stage failure at Tj_ 
Payload in orbit 1.0 * 1.0 AU, 1=2° 


RESCUE RESPONSE : 


Launch expendable OTV and baseline kick stage 

Use OTV for rendezvous and transfer AV 

Use kick stage for placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 

Time 

Distance 

C 3 or A V 

Inject to Earth escape 

T l + (N/2) y 

1.0 AU 

C 3 « 1.08 0 I = 2° 

Rendezvous and dock 

'vlO 1 * later 

1.0 

A V = 0.100 km/sec 

Transfer to 0.86 x 1.0 AU 

same 

1.0 

1.143 

Placement 

, v0.45 y later 

0.86 

1.187 

STAGE REQUIREMENTS: 




Stage Number 

aV (km/sec) 


Propellant (lbs) 

1 (Expendable OTV) 

4.656 

44687 (41645 used) 

2 (Baseline Kick) 

1.187 


8829 (7591 used) 

TOTAL INJECTED WEIGHT: 

51057 + 11943 = 
in 300 km Earth 

63000 lbs 
parking orbit 



30043 lbs to C 3 

!t 

■ 

O 

00 


TOTAL RESCUE MISSION TIME: 

0.45 years 
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Table 4-2 (cont'd.) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 


CASE NUMBER : 
FAILURE EVENT: 


1 (Alternative) 

OTV and kick stage failure at T[_ 
Payload in orbit 1.0 x 1.0 AU, 1 = 2° 


RESCUE RESPONSE : Launch two** baseline kick stages 

Establish best possible stable orbit, 
0.86 x 0.873 AU* 


RESCUE SCENARIO: 


Event 

Time 

Distance 

Cs or A V 

Inject to Earth escape 

T L + (N/2) y 

1.0 AU 

C 3 = 1.08 @ I = 2° 

Rendezvous and dock 

^10 d later 

1.0 

AV = 0.100 km/sec 

Transfer to 0.86 x 1.0 AU 

same 

1.0 

1.143 

Placement 

^0.45^ later 

0.86 

1.066* 

STAGE REQUIREMENTS: 

Stage Number 

AV (km/sec) 


Propellant (lbs) 

1 (Baseline Kick) 

0.875 


8829 (8829 used) 

2 (Baseline Kick) 

1.434 


8829 (8829 used) 

TOTAL INJECTED WEIGHT: 

2 x 11943 = 23886 

lbs to C 3 

11 

• 

O 

CO 


TOTAL RESCUE MISSION TIME : 0.45 years 


insufficient aV capability of 1.187 - 1.066 = 0.121 km/sec. 

**Second alternative using only single kick staqe is placement into circular 
orbit 0.915 x 0.915 AU. 
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Table 4-2 (cont’d.) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 


CASE NUMBER : 2A (Placement at T L + 2.344 y ) 

FAILURE EVENT : OTV and kick stage failure at T L 

Payload in orbit 0.93 * 1.0 AU, 1=2° 


RESCUE RESPONSE : Launch two* baseline kick stages for rendezvous, 

transfer and placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 

Time 

Distance 

C* or AV 

Inject to 0.859 x 1.0 AU 
orbi t 

Wait one revolution 

• 

+ 

—I 
1 — 

1.0 AU 

C 3 = 2.34 @ I = 2° 

Rendezvous and dock 

T l + 1.896 y 

1.0 

A V = 0.710 km/sec 

Transfer to 0.86 x 1.0 AU 

T l + 1 . 896 y 

1.0 

0.598 

Placement 

T l + 2 . 344 y 

0.86 

1.187 

STAGE REQUIREMENTS: 

Stage Number 

AV (km/ sec) 


Propellant (lbs) 

1 (Baseline Kick) 

1.083 


8829 (8829 used) 

2 (Baseline Kick) 

1.412 


8829 (8718 used) 


TOTAL INJECTED WEIGHT : 2 x 11943 = 23886 lbs to C 3 = 2.34 


TOTAL RESCUE MISSION TIME : 1.344 years 


♦Alternative using only single kick stage is placement into circular orbit 
0.90 x 0.90 AU, or elliptical orbit 0.86 x 0.941 AU. 
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Table 4-2 (cont'd. ) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 

CASE NUMBER : 2B (Placement at Tj_ + 4.240^ 

FAILURE EVENT : OTV and kick stage failure at Tl 

Payload in orbit 0.93 * 1.0 AU, I = 2° 

RESCUE RESPONSE : Launch two* baseline kick stages for rendezvous, 

transfer and placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 

Time 

Distance 

C* or AV 

Inject to 0.906 x 1.0 AU 
orbi t 

t l + 

1.0 AU 

C 3 = 1.61 01=2° 

Wait three revolutions 

- 


• 

Rendezvous and dock 

T L + 3.792 y 

1.0 

AV = 0.295 km/sec 

Transfer to 0.86 x 1.0 AU 

T L + 3.792 y 

1.0 

0.598 

Placement 

T L + 4 . 24 $ 

0.86 

1.187 

STAGE REQUIREMENTS: 




Stage Number 

AV (km/sec) 


Propellant (lbs) 

1 (Baseline Kick) 

0.893 


8829 (8397 used) 

2 (Baseline Kick) 

1.187 


8829 (7591 used) 

TOTAL INJECTED WEIGHT: 

2 x 11943 = 23886 

lbs to C 3 

= 1.61 


TOTAL RESCUE MISSION TIME : 3 . 240 years 


♦Alternative using only single kick stage is placement into circular orbit 
0.888 x 0.888 AU, or elliptical orbit 0.86 x 0.917 AU. 



Table 4-2 (cont'd.) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 


CASE NUMBER : 3A (Placement at T L + 2.242 y ) 

FAILURE EVENT : Kick stage failure at T L + 0.448 y 

Payload in orbit 0.86 x 1.0 AU, I = 2° 


RESCUE RESPONSE : Launch two* baseline kick stages for rendezvous 

and placement 0.86 AU circular 


RESCUE SCENARIO: 




Event 

Time 

Distance 

Cq or A V 

Inject to 0.715 x l .0 AU 
orbit 

t l + 1 . 6 * 

1.0 AU 

C 3 = 7.71 @1=2° 

Wait one revolution 




Rendezvous and dock 

T L + 1 . 794 y 

1.0 

AV = 1.460 km/sec 

Placement 

T l + 2 . 242 y 

0.86 

1.187 

STAGE REQUIREMENTS: 




Stage Number 

AV (km/ sec 

1 

Propellant (lbs) 

1 (Baseline Kick) 

1.460 


8829 (8808 used) 

2 (Baseline Kick) 

1.187 


6947 (6947 used) 

TOTAL INJECTED WEIGHT: 

11943 + 10061 = 

22004 lbs to 

C 3 = 7.71 


TOTAL RESCUE MISSION TIME : 1.242 years 


♦Alternative using only single kick stage is placement into elliptical orbit 
0.86 x 0.908 AU. 
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Table 4-2 (cont'd.) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 

CASE NUMBER : 3B (Placement at T L + 4.933 y ) 

FAILURE EVENT : Kick stage failure at T L + 0.448 y 

Payload in orbit 0.86 x 1.0 AU, 1=2° 

RESCUE RESPONSE : Launch one baseline kick stage for rendezvous 

and placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 

Time 

Distance 

Ca or AV 

Inject to 0.824 x 1.0 AU 
orbi t 

Wait four revolutions 

t l + 1 '° y 

1.0 AU 

C 3 = 3.20 @1=2° 

Rendezvous and dock 

T l + 4 . 484 y 

1.0 

AV = 0.429 km/sec 

Placement 

T l + 4.933 y 

0.86 

1.187 

STAGE REQUIREMENTS: 

Stage Number 

AV (km/sec) 


Propellant (lbs) 

1 (Baseline Kick) 

1.616 


8829 (8683 used) 


TOTAL INJECTED WEIGHT : 11943 lbs to C 3 = 3.2 


TOTAL RESCUE MISSION TIME : 3.933 years 


Table 4-2 (cont'd.) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 


CASE NUMBER : 
FAILURE EVENT: 


4A (Placement at T L + 2.989 y ) 

Kick stage failure at Tl + 0.448 y 
Payload in orbit 0.86 x 0.93 AU, I = 2° 


RESCUE RESPONSE: 


Launch one baseline kick stage for transfer, 
rendezvous and placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 

Time 

Distance 

C^ or AV 

Inject to 0.750 x 1.0 AU 

T, + 

1.0 AU 

C 3 = 5.88 01=2° 

orbit 

Transfer to 0.750 x 0.93 AU 

T, + 1 . 409 * 

0.750 

AV = 0.579 km/sec 

Wait 1.5 revolutions 

Rendezvous and dock 

T l + 2.565 y 

0.93 

1.191 

Placement 

T l + 2.989 y 

0.86 

0.622 

• 


STAGE REQUIREMENTS : 
Stage Number 
1 (Baseline Kick) 


AV (km/sec) 
2.392 


Propellant (lbs) 
8829 (8779 used) 


TOTAL INJECTED WEIGHT : 11943 lbs to C 3 = 5.88 

TOTAL RESCUE MISSION TIME: 1.989 years 
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Table 4-2 (cont'd.) 


RESCUE MISSION PROPULSION REQUIREMENTS - CASE STUDIES 

CASE NUMBER : . 4B (Placement at T L + 4.681 y ) 

FAILURE EVENT : Kick stage failure at T L + 0.448 y 

Payload in orbit 0.86 x 0.93 AU, 1=2° 

RESCUE RESPONSE : Launch one baseline kick stage for transfer, 

rendezvous and placement 0.86 AU circular 


RESCUE SCENARIO: 


Event 


Time 

Distance 

Cn or AV 

Inject to 0.805 x 1.0 

AU T. 

+ 1.0* 

1.0 AU 

C 3 = 3.76 @1=2° 

orbit 

L 




Transfer to 0.805 x 0. 

.93 AU T l 

+ 1 . 429 y 

0.805 

A V = 0.572 km/sec 

Wait 3.5 revolutions 





Rendezvous and dock 

t l 

+ 4.258 y 

0.93 

0.623 

Placement 

t l 

+ 4.681 y 

0.86 

0.622 

STAGE REQUIREMENTS: 





Stage Number 


A V (km/ sec) 


Propellant (lbs) 

1 (Baseline Kick) 


1.817 


8829 (7636 used) 

TOTAL INJECTED WEIGHT: 

11943 

lbs to C 3 = 

3.76 



TOTAL RESCUE MISSION TIME : 3.681 years 




Fig. 4-4 ORBITAL VARIATIONS FOR THE INITIAL ORBIT a - 0.90 AU, 







Earth Perihelion 



°* 70 0 r^ IE 05 2E 05 3E 05 4E 05 5E 05 6E 05 7E 05 8E 05 
, TIME (years) 

Fig. 4-6 ORBITAL VARIATIONS FOR THE INITIAL ORBIT a ■ 0.88 All, e = 0.022, 



Earth Perihelion 



Fig. 4-7 ORBITAL VARIATIONS FOR THE INITIAL ORBIT a = 0.89 AU, e « 0.0337 



for three slightly elliptical initial orbits: 0.86 * 0.88, 0.86 * 0.90, 

and 0.86 * 0.92 AU. These orbits also appear to be stable for periods up 
to 10 6 years or, at worst, marginally cross the Earth's orbit for brief 
intervals during this period. We would conclude, therefore, on the basis 
of these results obtained from secular perturbation theory, that alterna- 
tive options do exist for placement of failed payloads into off-nominal 
stable orbits. 

4.3 Automated Rendezvous/Docking Assessment 

The most critical aspect of rescue missions is undoubtedly the 
ability to rendezvous and dock with a payload in deep space. The rescue 
vehicle must be transfered to the near vicinity of the payload position, 
close enough to obtain radar acquisition of the target. It must then be 
very accurately guided during the terminal rendezvous phase through final 
closure with the target. Lastly, a secure dock must be accomplished in 
an automated mode. While these maneuvers are not yet commonplace in the 
space program, there does exist a background of experience and techno- 
logical development which allows one to project this as an engineering 
problem that can be solved. The Soviets have already demonstrated auto- 
mated rendezvous and docking in Earth-Moon space. The U.S. program in- 
cludes analysis and design experience related to backup options for 
manned space flight and, more recently, studies related to Mars sample 
return missions [1] and teleoperator control of Skylab's orbit decay 
problem [3]. 

Some of the more salient features and requirements of rescue opera- 
tions are outlined in Figure 4-8. The most important requirement is that 
the target vehicle should be able to play some cooperative role in the 
operations. Earth-based radio tracking of the target orbit for purposes 
of transfer guidance implies an operable communications link, at least 
in the transponder section of the payload's communication system. A 
transponder function is also needed to relay the rescue vehicle's radar 



Fig. 4-8 

AUTOMATED RENDEZVOUS AND DOCKING IN SOLAR ORBIT 
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signal during the terminal phase of rendezvous. Attitude control capa- 
bility through an active command/receive link with the rescue vehicle 
is certainly a vital requirement for success of the final closure and 
docking maneuvers. Target vehicle maneuverability other than attitude 
control is generally not necessary since the rescue vehicle is the active 
rendezvous partner. 

An estimate of target position knowledge error is shown in Figure 4-9 
assuming conventional DSN radio tracking. The error characteristic varies 
linearly with geocentric tracking distance and inversely with the sine of 
geocentric declination. A typical range of the RMS position uncertainty 
is 20 to 85 km at p = 0.5 AU and 40 to 170 km at p * 1.0 AU. Sensitivity 
to low values of declination can be reduced significantly by using quasi 
very long baseline interferometry (QVLBI) techniques in addition to con- 
ventional doppler tracking. These "new" data types involving simultane- 
ous (multi -station) doppler and range measurements are expected to be in 
common use for future planetary missions. Orbit transfer guidance accu- 
racy can be made to closely approach the orbit determination knowledge 
accuracy by employing several midcourse correction maneuvers to nullify 
AV execution error effects. Therefore, we might reasonably expect the 
3o error at terminal rendezvous initiation not to exceed 300 km. 

Radar acquisition of the target from several hundred kilometers 
should not present any problem to radar system design nor require exces- 
sive mass or power. Table 4-3 shows the characteristics of the rendez- 
vous radar proposed in the Martin-Marietta study of Mars sample return 
missions. The design, based on technology used in the Apollo rendezvous 
radar, is a unified S-band PM/CW system which serves a multipurpose func- 
tion for both rendezous and docking operations. It provides range, 
range-rate, and angle data from a maximum unambiguous range of 750 km 
down to a minimum docking range of 3 m. Due to rapidly degraded accuracy 
at very close range, range measurements are not used within 30 m of final 
docking; system performance relies instead on accurate range-rate and 






Table 4-3 


RENDEZVOUS RADAR CHARACTERISTICS (Martin-Marietta Ref.) 


System Parameters 

Frequency 

Radar Type .... 
Radar Mode .... 
Modulation . . . . 
Radar Power .... 
Maximum Range . . . 
Minimum Range . . . 
Radar Antenna . . . 
Angle Track Method 
Transponder Power . 
Transponder Antenna 
Coherence Ratio . . 


S-Band 

CW 

Automatic 

PM (819 kc Subcarrier, 4 minor tones) 
0.3 w (Solid State) 

750 km 
3 m 

Traveling Wave Array 
Phase Monopulse 
0.15 w (Solid State) 

Cassegrain 

220/239 


Error Summary 

Range Error (Bias) 3 m 

Range Error (Random) 

R < 65 km <3 m 

R = 750 km 750 m 

Range Rate Error (Bias) 5 cm/sec 

Range Rate Error (.Random) ...... 5 cm/sec 

Angle Error (Bias) 1.5 mrad 

Angle Error (Random) 

R < 10 km <0.05 mrad 

R = 750 km 3.2 mrad 
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angle data maintained through impact. Alternative systems such as a 
scanning laser radar have also been proposed for short-range applications. 
One important advantage would be to minimize the active RF interface with 
the target vehicle through the use of corner reflectors instead of a trans- 
ponder. Laser radar technology is currently in the breadboard stage and 
its future hardware development appears to be tractable. Additional sen- 
sor and control elements required on-board the rescue vehicle include 
celestial and inertial instrumentation for attitude control, dual -direc- 
tional thrusters for axial and lateral channel trajectory control, and a 
command sequencer/computer to implement the autonomous maneuver strategy. 

Although detailed engineering design and costing studies are clearly 
needed after the vehicle configurations are better defined, we would con- 
clude that rescue mission capability is technically feasible if a coopera- 
tive rendezvous mode can be assured . Such assurance implies a limited time 
interval between nominal payload launch and retrieval and, perhaps, some 
level of redundancy in payload vehicle. systems to enhance operational 
reliability. In the absence of cooperative rendezvous, new technology 
requirements are expected to be difficult to satisfy, particularly as 
regards capture mechanisms to implement automated docking with a massive, 
uncontrollable vehicle. One might consider an alternative risk reduction 
approach in such a circumstance— explosive destruction of the nuclear 
waste payload. Small remnant particles subject to erosion and ionization 
in the space environment can be strongly influenced by nongravitational 
forces such as solar radiation pressure and the solar wind. These forces 
tend to disperse the particle orbits and diminish the amount of material 
that would be intercepted by Earth. Those particles that do reenter 
Earth's atmosphere would most likely suffer complete burnup at high alti- 
tude into a submicron size oxide aerosol. The consequential risks of 
adverse health effects have been estimated to be quite small due to the 
extremely large dilution provided by the atmospheric volume of an entire 
hemisphere [2]. 
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